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foreword 


This study was initiated to determine the feasibility of predicting 
wake profiles behind buildings and natural obstacles using a scaled model 
in a wind tunnel. The wind tunnel approach is preferable because of 
economy of time and money, simplicity and convenience. This is the 
first report of a continuing program sponsored by the Fluid Dynamics 
Branch, Atmospheric Sciences Division of the Space Sciences Laboratory at 
the George C. Marshall Space Flight Center, National Aeronautics and 
Space Administration, Huntsville, Alabama. 

This research was conducted under the technical direction of Mr. Dennis 
W. Camp and Mrs. Margaret Alexander of the Space Sciences Laboratory at 
Marshall Space Flight Center. The support for this research was provided 
by Mr. John Enders of the Aeronautical Operating Systems Division, Office 
of Advanced Research and Technology, NASA Headquarters. 
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Coriolis parameter 
height of building or model 
von Karman constant 

constant in Couni han-Hunt-Jackson Theory calculated 
from Eq. 13 

subscript for model 

Momentum flow per unit width 

power law exponent 

subscript for prototype 

meteorological towers Nos. 1,2, 3, 4, 5, 6 

horizontal component of air velocity of any z 

horizontal component of air velocity at z = H upstream 
of building 

horizontal component of air velocity at z = 6i upstream 
of model 

horizontal component of air velocity far upstream of 
building (Tl) at z/H = 6.5 

geostrophic or free stream wind velocity 

velocity deficit based on upstream velocity 

friction velocity, 

RMS of fluctuation of horizontal component of air 
velocity at any z 

elevation 

roughness length 

model or prototype boundary layer thickness 
height of surface layer, i.e., logarithmic layer 
dimensionless velocity deficit defined by Eq. 12 
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n 


dimensionless height z defined by Eq. 11 


0 wind direction, 0 = 180° for wind from T1 to T6 

(wind from N 30° W) 

X geographic latitude 

V kinematic viscosity of the air 

p density of air 

T shear stress 

Tq surface shear stress 

0 ) earth's rotation rate 



CHAPTER 1 


INTRODUCTION 

The work described in this report is motivated by the need to under- 
stand the wind environment around air terminals. Ascent or descent of 
aircraft through the atmospheric boundary layer is accompanied by changes 
in lift associated with changes in wind speed with altitude. The accelera- 
tion produced by the brief action of unbalanced forces results in devia- 
tions from the original flight path for descending flight. The above 
effect is enhanced by induced flows produced by buildings or natural 
obstacles in the vicinity of airports. The effect of these surface 
obstacles on the aerospace environment around airports has been reviewed 
recently by Fichti, et al . [1]. Shear layers, or wakes, produced 
downwind of surface obstacles can prove hazardous to aircraft, especially 
those of the V/STOL type, because of the high rate of change of wind 
speed with altitude in the layer. Clearly research is needed to determine 
the locations of these regions of induced flows in the wakes of surface 
obstacles and their effects on aeronautical systems. 

Related theoretical and experimental research has been carried out in 
recent years. Most of the pertinent literature has been discussed 
in extensive literature surveys by Fichti, et al . [1] and by Frost 
[2]. Some of the work cited in these references appear to be useful in 
the present study. 

A particularly useful theoretical study has been reported by Counihan, 

et al. [3]. A theoretical framework is provided for correlating 
experimental wake measurements corresponding to a given building geometry 
and upstream profile characteristics. Models of long buildings immersed in 
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a boundary layer having a thickness of ten times the model height have 
been used by Woo, et al [4] to obtain measurements of wake velocity 
and turbulence. The data of Woo, et al [4] are expected to be useful 
in testing the predictions based on the theory of Counihan, et al [3]. 

Other experimental studies are also helpful in giving insight into 

the physics of the wake flow. Oka and Kostic [5] made detailed measurements 

of velocity and turbulence in the recirculation region behind a two- 

dimensional square rod placed on one wall of a channel of rectangular 

cross section. Mueller and Robertson [6] measured wake profiles behind 

two-dimensional obstacles of quarter-round cross section including both 

data in the recirculation region and far downstream. Good and Joubert [7] 

% 

studied the aerodynamic drag of bluff plates and obtained useful correla- 
tions between upstream conditions and plate height. These relations 
should be useful in predicting the effect of obstacle height on wake 
profiles. 

Most of the available reports of experimental research dealing 
with wake flow behind obstacles are based on work done with small 
scale laboratory equipment. Few field studies have been made, however, 
major studies have been reported by Frost and Shahabi [8], Frost, et. al 
[9] and Sacre [10]. In the work of Frost and Shahabi the wake of a 
simulated block building was studied under field conditions using ten 
instrumented wind towers located at Marshall Space Flight Center's 
Atmospheric Boundary Layer Facility (ABLF). The anemometers and vanes 
supported on the towers were used to measure mean horizontal wind profiles, 
wind direction and vertical wind speed. In addition to the three com- 
ponents of wind, the turbulence components were also measured. The 
tower locations are shown schematically in Fig. 1, where T1 , T2, T3, 
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T4, T5, and T6 are used to denoted the six principal towers used in the 
investigation. The distances given in Fig. 1 are in meters, and the 
dimensions of the simulated block building are indicated as height 3.2m, 
width 2.4m and length 26.8m. The investigation considered primarily 
winds in the direction from T1 to T6, i.e., winds from North 30° West. 

The last measuring station (T6) is located about 44 building heights 
downwind of the simulated building, i.e., x/H = 44 at tower T6. The 
instruments were located on the towers at 3, 6.2, 12 and 20.88m above 
the ground, i.e., at z/H equal to 0.94, 1.94, 3.75 and 6.5, respectively. 
Thus wind profiles up to 6.5 building heights above the ground are measured 
at several stations starting eight building heights upwind of the building 
(Tl) and extending to 44 building heights downwind of the building (T6). 

The boundaries of this wake study clearly encompass the principal regions 
of affected wind present in the flow. Too, these wind profiles are con- 
tained in the region traversed in the wind tunnel investigation of Woo, 
et al [4]. 

The wind tunnel investigation of Woo, et al [4] was carried out in the 
Meteorological Wind Tunnel (MWT) in the Fluid Dynamics and Diffusion 
Laboratory at Colorado State University (CSU) using several models of 
buildings, one of which was a 1/50-scale model of the 3.2 x 2.4 x 26.8m 
building depicted in Fig. 1 and used in the field study at MSFC-ABLF. 

The upstream profile was simulated by using the appropriate artificial 
roughness on the wall of the wind tunnel. Both mean and fluctuating 
velocities were measured in the wake of the model building. Profiles 
were measured at several downstream stations in the wake including four 
which correspond exactly to the locations of towers T4, S3, S4, and T5. 

Thus it is possible to use the wind tunnel data to compare with the actual 
field data taken in the atmospheric boundary layer. 
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The present work is a preliminary comparison of the MSFC-ABLF field 
data and the CSU-MWT wind tunnel data. The purpose of the comparison is 
to assess the accuracy of predicted wake profiles for the prototype 
building based on measurements made in the wake of the wind tunnel model, 
and to arrive at conclusions as to how disparities can be reduced or 
eliminated. Accurate prediction of building wakes for the purpose of 
mapping the wind environment of aircraft near terminals is most economically 
accomplished by wind tunnel modeling. However, evidence that wind tunnel 
wake data can be correlated with full-scale wake results is a necessary 
step in the establishment of the reliability of the technique. 
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CHAPTER 2 


COMPARISON OF UPSTREAM CONDITIONS 


This chapter is devoted to the characteristics of velocity profiles 
upstream of the building or model. Referring to Fig. 1 for the field 
tower arrangement, the upstream profile is that measured at the first 
tower. Tower T1 , assuming a wind direction from T1 to T6, i.e., the 
wind angle* 0 is 180°. The profile so determined can be characterized 
by the exponent n found in the power law 


This relation models the entire boundary layer, and the exponent n 
can be determined from velocity profiles measured either in the field 
or in the wind tunnel without knowing or 6. The exponent n is the 
slope of the straight line graph which represents U versus z on log-log 
paper. The surface or logarithmic layer has a velocity distribution which 
fi ts 


where u* and z^ are the velocity and length scales, respectively, for 
this layer. The friction velocity u* can also be determined from the 
velocity profile measured in the wind tunnel or the field without knowing 
the roughness length. The friction velocity is simply obtained from 
the slope of the plot of U versus In z. 


*where 0 is defined to be the direction of alignment of the towers and is 
not measured from north. The towers are aligned approximately 30° off 
(west) from north. 
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Table 1 consists of values of friction velocity for five tower 
locations and of upstream power law exponent determined by the least 
squares method. The results of these calculations for runs having a wind 
direction 0 of approximately 180°, i.e., north 30° west, are presented 
in Table 1. Some values are omitted because of nonlinearity or 
insufficiency of data. The origin of the x-coordinate is at the downwind 
face of the building; thus the towers T1 and T2 are located at x/H values 
of -8.65 and -1.84, respectively, and T4, T5 and T6 have x/H values of 
4.88, 16.44 and 43.94, respectively. Runs 8013-8038 were made without a 
building, and Runs 8407-8512 were made with the building in place. 

The upstream wind profile exponent observed in the field is clearly 
affected by the presence of the trees and bushes shown in Fig. 1 upwind 
of the first tower. Since tests were carried out over a period of 
three years and the size of the natural vegetation changed during this 
time, it is to be expected that the exponent n would change. The 1972 
values were 0.14-0.19, while the 1974 values varied from 0.24-0.31. 
Additionally, T1 was located 3.5m* (x/H = -2.93) from T2 in 1972 as 
compared with the 21.8m distance used in 1974 as shown in Fig. 1. 
Deviation of wind direction from 180° also produced variation of the 
exponent. The values of u* depend on the upstream surface roughness, 
the size of obstacles and the geostrophic wind velocity U^. Thus 
considerable variation of u* is observed at Tower T1 . 

A variety of related model testing was done in the CSU-MWT and is 
reported by Woo, et al [4]. The test condition and model 

*This distance applies to Runs 8013-8038 in Table 1. 
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Table 1. Friction Velocity (m/s) for a Wind Direction of 180 Degrees 


Run No. 

Exponent n 
at T1 

T1 

T2 

T4 

T5 

T6 

8D13 

0.14 

0.210 

0.323 

0.268 

0.283 

0.198 

8018 

0.14 

0.252 

0.302 

0.448 

0.299 

0.269 

8019 

0.19 

0.238 

0.298 

0.457 

0.251 

0.222 

8038 

— 

— 

-- 

0.242 

0.328 

0.225 

8407 

0.26 

0.428 

0.558 

0.583 

0.249 

0.222 

8408 

0.25 

0.388 

0.559 

0.566 

0.413 

0.223 

8501 

0.24 

0.510 

0.795 

0.382 

0.256 

0.407 

8502 

0.28 

0.635 

0.828 

0.762 

0.530 

0.431 

8503 

0.31 

0.723 

0.915 

0.784 

0.584 

0.485 

8504 

0.26 

0.552 

0.808 

0.809 

0.474 

0.423 

8512 

0.31 

0.575 

0.512 

0.604 

0.428 

0.368 



sizes used are presented in brief form in Table 2. The model tests of 
principal interest in the present comparison study are those for test 
condition No.' 2 for which n = 0.27, 6 = 0.61m, = 0.0061m and = 16m/s. 

A field run with an exponent close to this value is Run No. 8407 with 
n = 0.26 and u* = 0.428m/s. Equation (2) can be used to determine z^ 
for this field run. The resulting apparent roughness length is 0.16m, 
which is excessive. Profile data for other runs at T1 also yield 
excessive values; e.g., the calculated values of z^ are 0.226m, 0.254m 
and 0.591m for Run Nos. 8408, 8504 and 8512, respectively. This is an 
order of magnitude larger than the value determined by Frost, et al . [9], 
viz., z^ = 0.007m. The explanation for this follows from Equation (2) 
in the form 


= ovn kU 

0 6XP 


(3) 


coupled with the observation that an obstacle causes an initial increase 
followed by a decay of u*, as can be observed from values given in Table 1 
of this paper. Since the obstacles, i.e., trees, bushes, fences and 
ditches, upstream of Tower T1 initially raise the value of u*, and since 
u* decays very slowly, u* at Tower T1 is expected to be larger than its 
equilibrium value, and Equation (3) shows that the corresponding calculated 
value of Zq should also be higher. On the other hand. Equation (3) yields 
values of z^ of 0.0107m and 0.0116m at Tower 6 for Run Nos. 8504 and 8512, 
respectively, which are reasonable values of roughness length for the 
terrain used. 

It is clear that care must be exercised in modeling atmospheric flows 
by simply adjusting wind tunnel power law exponents. Although the velocity 
profiles for wind tunnel condition No. 2 and field Run No. 8407 have 
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Table 2. Geometric Condition for Modeling in CSU-MWT 


Test 

Condition 

6(m) 

ZQ(mm) 


Model 

1 

Model 

2 

Model 

3 

Model 

4 

H/6 


H/6 


H/6 


H/6 


1 

0.71 

4.3 

165 

0.044 

7.21 

0.0915 

15.1 

0.116 

19.2 

0.090 

14.9 

2 

0.61 

6.1 

100 

0.051 

5.08 

0.1066 

10.7 

0.135 

3.5 

0.105 

10.5 

3 

0.38 

0.11 

3455 

0.082 

281.8 

0.171 

591 

0.217 

751 

0.168 

582 

4 

1.37 

0.048 

28542 

0.023 

645.8 

0.047 

1354 

0.060 

1721 

0.047 

1333 



approximately the same characteristic exponent, i.e., n = 0.27, the 
upstream wind tunnel profile is an equilibrium profile deriving from the 
uniform surface roughness of the wind tunnel wall, but the field profile 
at Tower T1 is a non-equilibrium profile in which turbulence generated 
by upstream obstacles is still decaying in the flow direction. Thus 
there is reason to expect some dissimilarity in the wake profiles of 
the model and prototype at stations downstream of the building resulting 
from differences in upstream velocity profiles. Dissimilarity of model 
and prototype wake flow fields would also be expected if geometric or 
turbulence dissimilarities were present in the upstream flow. 
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CHAPTER 3 


CONDITIONS OF SIMILARITY 

Cermak [11] and Sundaram, et al [12] and Armitt and Couni han [13] 
have discussed the conditions of similarity for correct modeling of 
atmospheric surface layers. Among the several criteria for similarity 
of model and prototype flows is geometric similarity. For the present 
case this means equality in the ratios of upstream boundary layer 
thickness 6 to building height H or roughness length for model 
and prototype. Values of these length ratios used in the wind tunnel 
study of Woo, et al [4] are summarized in Table 2. Of interest in 
the present comparison are values for condition No. 2, model No. 4 
which were used in the wind tunnel investigation. These can be compared 
with corresponding values in the field study obtained by using the 
formula for estimating planetary boundary layer thickness given by 
Blackadar and Tennekes [14], viz., 

« = i; (4) 

If the value of u* at T6 of run No. 8407 from Table 1 (u* = 0.2218m/s) 
is taken as a near equilibrium value, and f is calculated from 

f = 2o) sin X (5) 

-4 -1 

as 0.8133 X 10 sec , then 6 is estimated to be 682m. These values 
yield prototype length ratios of H/6 = 0.0047 and H/z^ = 457 as compared with 
H/6 = 0.105 and H/z^ = 10.5 for model. Thus noticeable dissimilarities 
appear in relative heights of surface roughness elements and the obstacle 
with respect to overall boundary layer thickness, and it is clear that the 
building extends into a much smaller part of the atmospheric boundary layer 
than does the model in its boundary layer. However both obstacles are in 
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the logarithmic layer governed by Equation (2), since Plate [15] estimates 
the thickness of 6^ O'*" *^*^® logarithmic layer in a wind tunnel as 
0.156 > H = 0.1056, and Blackadar and Tennekes [14] give the atmospheric 
surface layer thickness, i.e., the patching height, as 0.046 > H = 0.00476. 
The ratios H/6|_ are unequal for model and prototype, but both upstream 
velocity profiles are logarithmic. The lack of geometric similarity is 
related to a corresponding dissimilarity in the turbulence structures 
of model and prototype flow fields. 

The root-mean-square of the longitudinal component of velocity fluc- 
tuation u' is shown in Fig. 2 for the wind tunnel (solid curve) and the 
field (triangular and square points) as a function of nondimensional 
height z/H. The data are seen to become more disparate with increasing 
z/H, although some agreement is shown for z/H < 2. The disparity at 
large z/H is certainly explained by the difference in relative depth of 
the boundary layers, since the ratio of (6/H) prototype to (6/H) model 
is 22.38. This effect is illustrated by consideration of the field data 
point at z/H =6.5, which would move to z/H = 0.29 if the 6/H ratios are 
taken into account. Thus the field data constitute only a very small 
segment of the u'/u* curve for the whole boundary layer. Field data 
in the range 2.0 < u'/u* < 2.5 agree with other field and wind tunnel 
data for the constant shear stress layers, e.g., see Armi.tt and Counihan 
[12]. However, complete agreement of wake velocity profiles for proto- 
type and model flows would not be expected in view of the turbulence 
dissimilarities at z/H > 2. 

In the above discussion the heights of the atmospheric boundary layer 

and the surface layer have been crudely estimated. The assumed similarity 

between model and prototype layers at levels above the surface layer does 

not really exist. In view of the crudeness of the model implied by Eq. 1, 
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Sundaram, et al [12] suggest that the height of the logarithmic 

layer 6^^ should be considered in lieu of 6 in the development 

of similarity parameters for the flow field around surface obstacles. 

Under this assumption the conditions for similarity become 

6 , 6 , 

^ (FT^P 


and, as previously used. 


(iL) = (M_) 


( 7 ) 


An additional condition to assure that the effects of molecular viscosity 
are negligible, i.e., to assure an aerodynamically rough tunnel wall, is 

u* 


(- 


V 




( 8 ) 


Combination of Equations (2), (6) and (8) leads to 


(< 5 , ) 


3v 


L^M - k U, 


M 


6 , 6 . 


(9) 


Equation (9) gives 6^^ = 0.32m for the minimum thickness of the logarithmic 
layer in the wind tunnel. Figure 2 shows that 2H = 0.128m for the 
actual wind tunnel flow. However, the required obtained from 
Equation (9) could be reduced by running the wind tunnel at a higher 
speed Uj^. In order to satisfy Equations (6) and (7), the model should 
have a height H of 1.5 cm, and the wind tunnel wall should have a rough- 
ness length of 0.0000355m. These values are based on 6 = 682m and 6|^ 

= 0.046 = 27.28m, and these values could conceivably be larger by a 
factor of three or four, thus increasing the model height H to the 
value actually used (6.4 cm). However, calculating 6 from Equation (4) 
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using the data for u* from Table 1 (T6), a factor of two, i.e., 6 1400m, 

is a more reasonable upper limit for the field runs listed. 

In general, the similarity conditions of Sundaram, et al [12] 
require a higher tunnel speed, a smaller model size and a smoother 
tunnel wall. Since the wind tunnel conditions actually used by Woo, 
et al [4] to model the field conditions studied by Frost and Shahabi 
[8] are not strictly similar, some disparity is to be expected in the 
wake profiles. On the other hand, at least qualitative similarity of 
wake profiles should be expected owing to the approximately equal and 
constant (with height) upstream values of u'/u* in the logarithmic 
layer and equal power law exponents n for upstream velocity profiles. 
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CHAPTER 4 


COMPARISON OF WAKE VELOCITY PROFILES 

The effect of the building or model on the upstream velocity profile 
is studied in this chapter by comparison of field and wind tunnel data 
presented graphically in non-dimensional form. Figure 3 shows the wind 
tunnel velocity profile for Test Condition 2 upstream of Model 4 compared 
with field profiles at tower T1 for Runs 8407, 8504 and 8512. In this 
figure the reference velocity, used is that at tower T1 , level 4, 
i.e., at z/H = 6.5. Data for Run 8407 agree almost perfectly at all 
levels except z/H = 1.94. Data for Run 8504 make up a slightly fuller 
profile but are in excellent agreement, and those for Run 8512 indicate 
a slightly straighter profile. 

Profiles at tower T4 are compared in Fig. 4 for the two field runs 
8407 and 8504. The profile of Run 8504 appears to be slightly shifted 
to the right, which may indicate some unaccountable external effect 
which has increased the layer's momentum at every level over that observed 
in Run 8407. The possibility that this effect is produced by a change 
in height of natural obstacles (hedges and trees) upwind of the 
simulated building associated with the eight-month difference in time 
between Run 8407 and Run 8504 must be eliminated because of the very good 
agreement at level 4 of the data of Run 8512 with the earlier data. 

The enhancement of momentum at level 4 is, however, also observed 
for Run 8407 far downwind in the wake, i.e., at tower T6. This may be 
inferred from a comparison of data for Runs 8407 and 8512 in Figs. 4 and 6 
Unfortunately, good data at level 4 tower T6, is not available, but 
it is clear that data for 8407 follows 8512, and thus that significant 
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momentum enhancement is evident at tower T6 for both of these runs, as 
well as Run 8504. 

A rough estimation of momentum flow per unit width MOM is obtained 
by evaluating the integral 

1 - 20.88 5 

MOM = U^^dz (10) 

0 

3 

for Run 8504. This integral is evaluated at 998m /s at T6 compared 

3 

to 629m /$ at T1 . This is an increase of 59% over the upstream momentum 
for the layer defined by 0 < z < 20.88m. However, the wind tunnel data 

3 

applied to the corresponding model layer yields 47m /s at T6 compared 

3 

with 57.3m /s at T1 , i.e., a momentum decrease of 18%. Since a momentum 
decrease is expected, the validity of the field data is suspect. However, 
in view of the conclusions reached in Chapter 2 of this report, the 
upstream field profile is not an equilibrium profile, whereas the wind 
tunnel profile is a fully developed profile, i.e., the field profile is 
in the process of receiving horizontal momentum from the layer above it 
at the time it passes over the simulated building (and loses momentum) 
and the process of momentum addition continues in the wake of the building 
with a condition of near equilibrium being achieved at station T6. Thus 
the fullness of the profile at T6 would be present at T1 also if equilibrium 
had obtained there. Data from other runs were used to evaluate the momentum 
integral of Equation (10), and the results are presented in Table 3. These 
data provide convincing evidence that the wind profiles at T1 are not 
equilibrium profiles. 

The possibility that the apparent increase in momentum is a momentum 
decrease at T1 associated with the presence of the building 8.65 building 
heights downstream of T1 . Based on the measurements of Rider [16] it is 
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estimated that this effect could indicate an anomalous increase of as 
much as five percent. However, values shown in Table 3 indicate much 
larger increases. 

Table 3. Momentum Integral Evaluation 



q 

(Units are m /s) 


Run No. 

MOM at T1 

MOM at T6 

8501 

554 

916 

8502 

660 

919 

8503 

650 

1007 

8512 

421 

617 

8407 

334 

428 

8408 

288 

410 


Figures 3 through 6 show what is though to be a continuing adjustment 
of field wind profiles to local surface roughness combined with a recovery 
from the retardation associated with the passage of the air over the 
building. Since the wind tunnel profile is merely returning to equilibrium 
and its value is greater than the corresponding field value, the wind 
tunnel profile shown in Fig. 6 does not match that measured in the field 
for a much smaller relative roughness. The disparity between model and 
prototype wake profiles is illustrated qualitatively in Fig. 7 where Equa- 
tion(l) is used to determine equilibrium profiles using n = 0.27 for the 
wind tunnel and n = 0.13 for the field. The latter exponent is reasonable 
for Zq = 0.7 cm (e.g., see Plate [15], p. 41). Thus qualitative prediction 
of the downwind wake profiles (x/H = 40) is possible a priori, given only 
the values of z^ for the prototype and model. 
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The relative effect of the obstacle on the flow may also be observed 
in Figs. 3 through 6. A very significant retardation is noted in the 
wind tunnel wake, but almost no retardation is observed in the field 
profiles. Accompanying the greater retardation is a 75% higher velocity 
gradient at the lower levels which should lead to greater turbulence 
levels in the wind tunnel wake. These observations can be explained 
by the greater value of H/6 in the wind tunnel flow. Figure 8, which is 
based on the correlations of Good and Joubert [7], shows qualitatively 
that the effect of increasing H/6 is to increase Cp for the obstacle. 

The graphs indicate that the use of a model with H/6 greater than the 
prototype and with H/z^ less than the prototype would produce a larger 
Cjj, and hence a greater momentum deficit and turbulence excess would 
be expected. 
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CHAPTER 5 


COMPARISON OF WAKE TURBULENCE PROFILES 

Figure 2 was used to show the comparison of upstream turbulence 
distribution. Figures 9 through 15 extend the comparison to wake pro- 
files. Figures 9 through 12 show the RMS value of longitudinal tur- 
bulence u' nondimensionalized with reference mean velocity U^. The up- 
stream profiles for model and prototype shown in Fig. 9 do not intersect 
at any point. This is a different behavior than previously exhibited 
in Fig. 2 for u'/u* profiles using the same turbulence data. On the 
other hand, the field data plotted in Fig. 9 are in close agreement. 

The same apparent independence of field and wind tunnel data is ob- 
served at the other stations shown in Figs. 10 through 12. The near 
equilibrium profiles of Fig. 12 are very close to the starting profiles 
of Fig. 9. 

The response of the flow to the obstacle shown in Fig. 10 is stronger 
for the model than for the prototype, e.g., u' increases 70% for the model 
and 15% for the prototype at level 2. This effect agrees with the difference 
in velocity gradient observed in Fig. 4 for the same station, i.e., the 
wind tunnel velocity gradient is 75% higher than the field gradient 
in the lower part of the layer, which implies a correspondingly higher 
relative turbulence production for the wind tunnel layer just downstream 
of the obstacle. Another point of agreement is the previously noted 
difference in momentum flow, i.e., the greater increase in velocity grad- 
ient observed in Fig. 4 is associated with a greater momentum loss, which 
implies a greater relative drag force on the obstacle. A greater relative 
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Fig. 12. Turbulence Profiles Far Downstream 
at x/H = 43.94 
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drag force is expected for the greater value of H/6 present in the wind 
tunnel as was shown in Fig. 8. 

Turbulence wake profiles, non-dimensional ized with friction velocity, 
are shown in Figs. 13 through 15. More scatter of field data points 
is observed in these graphs, but a better mingling of wind tunnel and 
field points is noted. This is particularly true at the farthest down- 
stream station as shown in Fig. 15. Here wind tunnel data for levels 
one through three agree almost perfectly with the field data. 
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CHAPTER 6 


COMPARISON WITH WAKE THEORY 

An alternative to direct comparison of velocity profiles is a 
comparison of model and prototype profiles with those obtained from 
wake theory. One way this can be done is to calculate an approximate 
value of the parameters introduced by Counihan, et.al. [3], which 
are called r, and n herein, and are defined as 


- L ^KXx 
^ H ^H ^ 


n+1 


and 


^ - U„H 


with the constant K calculated from 

U. 


K 


-1 


2k u* 


( 11 ) 


( 12 ) 


(13) 


The velocity deficit AU is defined by Frost and Shahabi [8] as the loss 
in momentum per unit mass at a given height, i.e., U at a downstream 
station minus U upstream of the building at the same distance z above 
the ground. The wind tunnel profiles shown in Fig. 16 are based on the 
foregoing definition of AU. However, the field profile, indicated with 
open square points, is based on defining AU as U at T5 minus U at the 
same level of T6. This is done since the profile at T1 is apparently 
not fully developed, so that its use would result in AU > 0 in the wake. 
Although the alternate definition of AU is necessary and useful, it is 
impossible to draw conclusions regarding the comparison of the field 
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profile with theory or with the wind tunnel profiles in the deficit 
form. 

The wind tunnel measurements shown in Fig. 16 are qualitatively 
similar to those presented by Counihan, et.al. [3]. A theoretical 
curve, calculated by the methods of Counihan, et.al. [3] is also 
included in Fig. 16 for comparison with wind tunnel wake profiles. 
Approximate agreement of theory and measurement is observed for x/H = 20 
and below. A large disparity is noted at x/H = 40. Below n - 1.5 
the correlation is extremely good up to station x/H = 25, indicating 
a self-preserving s-n relation within these limits. Above n 1.5 
the ?-ri profiles tend to spread somewhat, so that the profile depends on 
x/H.' Attempts are being made presently to modify the calculated 
parameters to produce a better correlation. 

Attempts to correlate data for the purpose of predicting the wake 
velocity and turbulence profiles have been investigated. Some results 
are presented in the appendix. Correlations of this sort should be 
very useful tools in predicting full scale wake behavior. 
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CHAPTER 7 


CONCLUSIONS 

Several conclusions may be drawn from a preliminary comparison 
of field data taken from the MSFC-ABLF and wind tunnel data obtained 
from the CSU-MWT. Similarity conditions require a higher tunnel speed, 
a smaller model size and a smoother wind tunnel wall. The observed 
disparity in wake wind profiles stems from differences in 6/z^ which 
are associated with differences in power law exponents. The percent 
turbulence increase generated by the model (70%) is larger than the 
corresponding u' increase (15%) developed by the prototype; this 
difference may be accounted for by the greater value of H/6 present 
in the wind tunnel. Field turbulence profiles, u'/u* as a function 
of z/H, are roughly predictable from model data at x/H = 44 provided 
z/H < 3. Wake theory agrees approximately with wind tunnel data for 
x/H £ 20. A self preserving wind tunnel wake, indicated by a good 
correlation of data in a ?-n plot exists for n <. 1.5 and x/H < 25. 

Further work is needed to achieve adequate data comparison. 
Suggested studies would include the reverse flow case (0 = 0°), v' 
and w' turbulence components, velocity deficit decay, turbulence decay 
and spectra, three dimensional effects, recirculation zone data, and the 
perturbation of non-equilibrium flow. 
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APPENDIX 


The internal boundary layer concept has been applied to wake flow. 

Fig. A1 shows upper and lower knees in wake velocity profiles. These 

knees mark the edge of regions of influence; the upper knee marks 

the outer edge of the internal boundary layer and the lower knee 

the edge of the sublayer. Fig. A2 shows a correlation of these data; 

1/2 

the layers grow according to x . 

A jet type correlation is used for velocity profiles near the 

building, and it is seen from Fig. A3 that the data are correlated 

by a single curve. Here ? is the usual jet parameter a z'/x. A similar 

correlation of turbulence data was attempted in Fig. A4. For the latter 

correlation to be useful a way of predicting u' must be found. 

max 

The decay of the velocity deficit is shown in Fig.A5. Presently 
a model is being developed to explain the change in slope of the decay 
curve. 
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Fig. A2. Growth of Internal Boundary Layers 
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Fig. A3. Velocity Correlation in Near Wake 
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Fig. A4. Turbulence Correlation in Near Wake 
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